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Hard photons emitted from energetic heavy ion collisions are very interesting since they do not 
experience nuclear interaction, and therefore they are useful to explore properties of nuclear matter. 
We investigated hard photon production and its properties in intermediate energy heavy-ion colli- 
sions with the help of the Blotzmann-Uehling-Ulenbeck model. Two components of hard photons 
are discussed: direct and thermal. The positive directed flow parameter and negative elliptic flow 
parameter of direct photons are demonstrated and they are anti-correlated to the flows of free pro- 
tons. The dependencies of hard photon production and anisotropic parameters on impact parameter, 
beam energy, nuclear equation of state and symmetry energy are also discussed. Furthermore, we 
investigated the two-photon momentum correlation function from which the space-time structure 
information of the photon source could be extracted as well as the two-photon azimuthal correlation 
which could provide another good method to determine the elliptic flow parameter V2 of direct hard 
photons. 

PACS numbers: 25.75.Ld, 24.10.-i, 21.60.Ka 
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I. INTRODUCTION 

One of the main goals of intermediate energy heavy 
ion collision (HIC) is to study properties of nuclear mat- 
ter, especially to determine the nuclear Equation-of-Statc 
(EOS). HIC provides a unique means to compress nu- 
clear matter to hot and dense phase within a labora- 
tory environment. The pressures that result from the 
high densities achieved during such collisions strongly 
influence on the motion of ejected matter and provide 
the sensitivity to the EOS. In comparison with the 
conventional hadronic probes, photons interacting only 
weakly through the electromagnetic force with the nu- 
clear medium are not subjected to distortions by the fi- 
nal state (neither Coulomb nor strong) interactions and 
therefore photon delivers an undistorted picture of the 
emitting source [l|-|jl- 

Since the last two decades, many model calculations 
and experimental facts [l|-|9[ have indicated that in inter- 
mediate energy heavy-ion collisions hard photons defined 
as photons with energies above the giant dipole reso- 
nance domain, above 30 MeV in this paper, mainly origi- 
nate from incoherent proton-neutron bremsstrahlung col- 
lisions, namely p-\-n — >■ p-\-n + ^. A nice review paper on 
hard photon was available [3|. These hard photons are 
emitted from two distinct sources. The first and dom- 
inant component denoted as direct hard photon (called 
as "direct photon" for short) is associated with the first- 
chance proton-neutron collisions in the initial phase of 
the heavy-ion reaction. The second one originates from 
the secondary proton-neutron collisions in the later stage 
of the reactions when the di-nuclcar system tends to be 
thermalized, accordingly called as thermal hard photon 
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(called as "thermal photon" for short). Because of their 
distinct emission sources, direct photons and thermal 
photons can deliver thermodynamic and dynamical in- 
formation of hot and dense nuclear matter formed during 
the various stages of the heavy-ion collisions. 

In relativistic heavy ion collisions, photons are also of 
very interesting since they can be served as one of the 
potential signals of quark-gluon plasma (QGP) forma- 
tion, eg., see Ref. |10l - KL3| . A hot QGP radiates a large 
amount of thermal photons, which dominate the spectra 
at small transverse momenta, whereas hard processes in 
nucleon-nucleon scatterings produce large momenta pho- 
tons. Therefore photon enhancement at low transverse 
momenta could be seen as a QGP signal, which has been 
observed at the BNL-Rclativistic Heavy- Ion Collider [l^l . 
Similar to the feature of photon production in intermedi- 
ate energy HIC, photons emitted from the interior of the 
hot matter no longer interact with the hadronic medium, 
in contrast to hadronic observables. 



The paper is organized as follows. In Sec. II, the 
simulation tool which we used is briefly introduced and 
the calculation method of photon production is outlined. 
Sec. HI describes the classification of hard photons and 
definition of anisotropic flow, and presents the results of 
the azimuthal asymmetry for direct photons and free pro- 
tons. In Sec. IV, we discuss the different variables (im- 
pact parameter, beam energy and EOS) dependences of 
the hard photon production and/or anisotropic flow pa- 
rameters. Sec. V gives the results and discussions on two- 
photon correlation functions, namely momentum corre- 
lation function, which is also called Hanbury-Brown and 
Twiss (HBT) correlation function, as well as azimuthal 
correlation function. Finally the summary is given in Sec. 
VI. 



II. BRIEF INTRODUCTION TO THE MODEL 
AND PHOTON PRODUCTION 

A. The Blotzmann-Uehling-Ulenbeck Equation 

The transport model is very useful for treating heavy- 
ion collision dynamics and obtaining important informa- 
tion of nuclear matter. In intermediate energy heavy- 
ion collisions, the Blotzmann-Uehling-Ulenbeck (BUU) 
model is an extensively useful tool [ij . The BUU equa- 
tion takes both Pauli blocking and the mean field into 
consideration, reads 



energy conservation. The double differential probability 
is given by 
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Here R is the radius of the sphere, ac is the fine structure 
constant, Pi and /3/ arc the initial and final velocity of 
the proton in the proton-neutron center of mass system, 
and 6 J is the angle between incident proton direction and 
photon emitting direction. More details for the model 
can be found in Ref. [l3, H • 
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It is solved with the method of Bcrtsch and Das Gupta 
[Tsj . In Eq.([T]), ^1^ and Vu are in-medium nucleon- 
nucleon cross section and relative velocity for the collid- 
ing nucleons, respectively, and U is the mean field poten- 
tial including the isospin-dependent term: 



U{p, r,) =. fl(^) + b{^r + Csym ^^" ^"^ T., (2) 

Po Po Po 

where po is the normal nuclear matter density; p. p„, 
and pp are the nucleon, neutron and proton densities, 
respectively; r^ equals 1 or -1 for neutrons and protons, 
respectively; The coefficients a, b and a are parameters 
for nuclear equation of state. Three sets of mean field 
parameters are used, namely the soft EOS with the com- 
pressibility K of 200 MeV (a= -356 McV, b = 303 McV, 
a = 7/6), and the semi-soft EOS with K of 235 MeV 
(a = -218 MeV, b = 164 MeV, a = 4/3), and the hard 
EOS with K of 380 MeV (a = -124 MeV, b = 70.5 MeV, 
CT = 2). Csyrn is the symmetry energy strength due to 
the density difference of neutrons and protons in nuclear 
medium, which is important for asymmetry nuclear mat- 
ter (here Csym = 32 MeV is used), but it is trivial for the 
symmetric system studied in the present work. 



B. Production cross sections of bremsstrahlung 
hard photon 

The BUU model was shown to be very successful in de- 
scribing the bulk properties of the reaction and nucleon 
emission in intermediate-energy heavy-ion collisions. In 
addition, the proton-neutron bremsstrahlung photon can 
be simulated as well in the model. For determining the 
elementary double-differential hard photon production 
cross sections on the basis of individual proton-neutron 
bremsstrahlung, the hard-sphere collision was adopted 
from Ref. [19, and modified as in Ref. 171 to allow for 



HI. AZIMUTHAL ASYMMETRY OF HARD 
PHOTONS 

A. Definitions of direct photons and thermal 
photons 

In the present work, we simulate the reaction of '^''Ca 
-I- "'"Ca collisions in most cases. Sometimes "'^Ca + "'^Ca 
and Kr + Ni are also simulated. 

Fig. [TJa) and (b) show the time evolutions of the pro- 
duction rate of bremsstrahlung hard photons and of sys- 
tem densities, including both maximum density (closed 
circles) and average density (open circles), respectively, 
for ^°Ca + ^°Ca colhsions at 60A McV in the ccntrality 
of 40%-60%. It is found that the hard-photon produc- 
tion rate is sensitive to the density oscillations during the 
whole reaction evolution. With the increase in density 
when the collision system is in the compression stage, the 
system produces more hard photons. In contrast, when 
the system starts to expand, the hard photon production 
rate decreases. Actually, the density oscillation of the 
colliding heavy-ion systems can be observed in the ex- 
periments via hard-photon interferometry measurements 
\M,M uM ^s well as dynamical dipolc 7 radiation [2fl- 
|22|. Apparently, hard photons are mostly produced at 
the early stage of the reaction. Thereafter we call these 
photons, emitted before the time of the first maximum 
expansion of the system (i ^ 65 fm/c in Ca -I- Ca at 60^ 
McV), as direct hard photons (on the left side of the blue 
dashed line in Fig. [ija)). It is also coincident with the 
definition in the Sec. I of direct hard photons. We call 
the residual hard photons producing in the later stage as 
thermal hard photons (on the right side of blue dashed 
line in Fig.[T|). In this way, we can identify the producing 
photon as cither direct or thermal by the emission time 
of photons in the present simulation. We notice that the 
production rate of thermal photon tends to decrease for 
later time, after 300 fm/c for Ca + Ca at 60^ McV, 
which can be understood by the decrease of n-p colli- 
sions in an expanding system. However, the hard photon 
yields in the later stage after 300 fm/c is relative small 
fraction of the total thermal hard photon emission. Also, 
in the following calculations, we mainly focus on direct 
hard photons which domintc the yield and do not depend 
on the time evolution for the later stage. 
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FIG. 1: (Color online) (a) Time evolution of hard photon produc- 
tion rate for the reaction '"'Ca + '"'Ca at 60^1 MeV in 40%-60% 
centrality. The EOS parameters with the compressibility K of 235 
MeV are used; (b) Time evolution of reduced maximum density 
Pmax/ po (closed circlcs) and reduced average density {p)/po (open 
circles) of the whole reaction system in the same reaction. The 
blue dashed line represents the time when till the first expansion 
stage, and in the panel (a) it separates direct hard photons (on the 
left side) and thermal hard photons (on the right side). 



For an example, Fig. [2] demonstrates that the direct 
and thermal hard photons for 60 A MeV -^OCa + ''°Ca 
collisions with a compressibility K of 235 exhibit differ- 
ent spcctuml shapes: thermal photons give rise to a softer 
energy spectral than direct ones. In the nucleon-nucleon 
cente-of-mass system, the photon spectrum can be de- 
scribed by the function: 



da 
IK 



= Iqc 



(4) 



where /q is a normalization constant and Eq is the slope 
parameter, which reflects the apparent source tempera- 
ture of the photon emission source. The slope parameter 
Eq depends on the bombarding energy and on the aver- 
age intrinsic momentum of the participant nucleons. As 
less energy is available on average in second-chance n- 
p collisions than in first-chance collisions because most 
of the projectile kinetic energy is damped, the thermal- 
photon spectrum becomes much softer. If one adds the 
two sources together in the photon spectruin, one obtains 
an empirical photon spectrum of the form [ll, H, Isl : 
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with the constant K^.th defined by: 
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Id,th represents the intensity of each photon source, d 
stands for the direct photon and th for stands for the 
thermal photon. However, we notice that the change 
in slope of the photon yield could be also affected by 
the i/E^ factor which enters the elementary np — npj 
bremsstrahlung probability |2^. In this case, Eq. [S] 
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FIG. 2: (Color online) Kinetic energy spectrum of direct photons 
and thermal photons for inclusive events of *''Ca -|- *'^Ca at 60 A 
MeV with the compressibility K of 235 MeV. Black dots and blue 
dots represent the direct photons and thermal photons, respec- 
tively. Lines represent the exponential fits with th Eq. \5\ 



should be changed into 
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In the next paragraph, we will compare the difference for 
the fitted slope with Eq. [5] and Eq. [71 

As shown in the inset of the figure, the slope parame- 
tesr of direct photons and thermal photons are 16.7 and 
8.1 MeV, respectively, for inclusive events of "'"Ca -|- '*"Ca 
collisions at 60 A MeV with an incompressibility K of 235 
MeV if we fit the spectrum with the Eq. \5\ However, 
Eq. [7] will give a larger slope value, namely 22.5 and 9.8 
MeV, respectively, which is 35% and 20% higher than 
the fits without the 1/E^ factor (not plotted in the fig- 
ure). With two different fit formulas, we put the impact 
parameter dependence of the slope parameter together 
in Fig. 131 A remarkable difference is seen between the 
slope values with and without the 1/E^ factor, i.e. Eq. [71 
gives a higher apparent photon temperature than Eq. [SI 
In some previous experimental analysis [ll S ImI, Eq- [SI 
was mostly used to extract the apparent slope, but this 
is probably not correct since a factor 1/E-y in front of 
the exponential form can be derived from their elemen- 
tary cross section [see Eq. [3]. On the other hand, the 
slope of hard photons displays a slight decreasing be- 
havior with increasing impact parameter, indicating that 
the emission source becomes cooler in peripheral colli- 
sions. Moreover, direct photons are hotter than thermal 
photons as we expected. 



B. Definition of anisotropic flows 

It is well known that collective flow is an important 
observable in heavy ion collisions and it can carry some 
essential information on the nuclear matter, such as the 
nuclear equation of state |24J34| . Anisotropic flows are 
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FIG. 3: (Color online) The slope parameters of direct and thermal 
hard photons as a function of impact parameter for ^'^Ca+'^'^Ca at 
60A MeV with the incompressibility K of 235 MeV. Open symbols 
represent the results from the Eq. [7] and solid symbols are from 
Eq.m 



defined as different nth harmonic coefiicients Vn of an 
azimuthal Fourier expansion of the particle invariant dis- 
tribution lii 



dN 
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oo 
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where (f) is the azimuthal angle between the transverse 
momentum of the particle and the reaction plane. Note 
that in the coordinate system the z-axis is along the beam 
axis, and the impact parameter axis is labeled as the x- 
axis. 

The first harmonic coefficient vi represents the di- 
rected flow, 



Vl 



{cos(f>) 
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Px + Py is ^^^ transverse momentum. The 



where pt - 

second harmonic coefficient V2 represents the elliptic flow 
which characterizes the eccentricity of the particle distri- 
bution in momentum space. 



V2 - (cos(2</))) = ( 
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C. Directed and elliptic floAvs of direct photons and 
free protons 

In relativistic heavy-ion collisions directed and ellip- 
tic flows of hard photons have been recently reported 
in the experiments and through theoretical calculations 
[35l - l39j . demonstrating a very useful tool for exploring 
the properties of hot dense matter. However, there are 
no experimental data available so far on the flow of hard 
photons in intermediate energy heavy ion collisions. In 
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FIG. 4: The time evolution of elliptic flow (1)2) of hard photons 
(a) and directed flow {v\) (b) for 60A MeV Ca -|- Ca collisions at 
centrality of 40-60%. 



this context, here we calculate directed transverse or el- 
liptic flows in intermediate energy heavy ion collisions. 
In addition, given that hard photons mostly originate 
from brcmsstrahlung by individual proton-neutron colli- 
sion and that free nucleons are also emitted from nucleon- 
nucleon collisions, therefore it will be interesting to com- 
pare flows of protons and photons. 

Fig. |4] shows the time evolution of the average values 
of directed flow and elliptic flow of the photons. Con- 
sidering the nearly symmetric behavior for directed flow 
versus rapidity, here we calculate the average vi over only 
the positive rapidity range, which can be taken as a mea- 
sure of the directed flow. 

From Fig. 4, the directed flow of photons show rapid 
rising with positive values during the compression stage 
and later on it decreases to even negative value. Af- 
terwards, the directed flow remains negative since the 
system is never compressed. For elliptic flow, the behav- 
ior shows contrary trend to the vi and later on it shows 
oscillation for thermal photon emission. The times cor- 
responding to the peak or valley values of flows roughly 
keep synchronized with the compression or expansion os- 
cillation on the system evolution as shown in Fig. 1(b). 

From the above calculations, we learn that direct pho- 
tons are preferentially emitted in the out-of-plane (nega- 
tive ^2) direction and thermal photons arc emitted from 
a thermalizing system which makes their emission more 
isotropic (i.e. the oscillated elliptic flow) than the direct 
ones produced in the pre-equilibrium stage. In addition, 
thermal photons contribute less than 30% of the total 
yield of hard photon in the present model. Therefor we 
will focus on direct hard photons to discuss the flow re- 
sults. 

The quantity of directed transverse fiow at mid- 
rapidity can be also defined by the slope: F = 



d(v) 



, where {y)c m is the rapidity of parti- 
(•y)c.™=o 

cles in the center of mass and (px) is the mean in-plane 
transverse momentum of photons or protons in a given 
rapidity region. In Fig. [SJa) and (b), we show (p^) plot- 
ted versus the cm. rapidity Tjcm. for direct hard photons 
(a) as well as (px) plotted versus the reduced cm. rapid- 
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FIG. 5: (Color online) (a) Average in-plane transverse momentum 
of direct hard photons as a function of cm. rapidity for 60^1 MeV 
Ca + Ca collisions in the semi-central centrality of 40 — 60%. The 
dashed line segment is a fit over the mid-rapidity region —0.5 < 
Vcm. < 0.5. (b) Same as the left panel but for free protons. The 
dashed line segment is a fit over the mid-rapidity region —0.5 < 

{y/ybeam)<^.m,. < 0.5. 
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FIG. 6: (Color online) (a) and (b) is the azimuthal distribution of 
direct photons and of free protons, respectively, for 60^4 MeV Ca 
-I- Ca collisions in the semi-central centrality of 40 — 60%. Both of 
them are fitted to the 4-th order Fourier expansion. 



ity {y/ybearn)c.m. for frcc protons (b) for 60A MeV Ca + 
Ca coUisfons in the semi-central centrality of 40 — 60%. 
The EOS parameters with the compressibility K of 235 
MeV are used. The errors shown are only statistical. For 
emitted proton (free proton), we identify it in the BUU 
calculation as its local density p < l/8po. We take the 
values of flows when the system has been already in the 
freeze-out time at 180 fm/c. A good linearity is seen in 
the mid-rapidity region (—0.5, 0.5) and the slope of a lin- 
ear fit can be defined as the magnitude of directed trans- 
verse flow. The extracted value for the directed trans- 
verse flow of direct hard photons is about -1-3.9 MeV/c, 
and that of free protons is about —23.7 MeV/c. Thus 
direct hard photons do exist the directed transverse flow 
even though the absolute value is smaller than the pro- 
ton's flow, and the sign of its flow is just opposite to that 
of free protons. 

In order to extract the value of elliptic flow and reduce 
the error of flts, we fit the azimuthal distribution to the 
4i/i order Fourier expansion. Fig |6] shows an example for 
such fits with the elliptic flow parameter V2= -2% for 
photons and V2~ 5.5 % for protons. 

Fig. [7] shows the differential elliptic flow of direct hard 



FIG. 7: Elliptic flows (1)2) for direct hard photons (a) and free 
protons (b) as a function of transverse momentum (pt) for 60A 
MeV Ca -|- Ca collisions in the semi-central centrality of 40 — 60%. 



photons (a) and of free protons (b) as a function of trans- 
verse momentum pT for *'^Ca -I- *°Ca at 60^ MeV. Sim- 
ilar to the directed transverse flow, the values of elliptic 
flow of direct hard photons and of free protons also have 
an opposite signs at this reaction energy, i.e., reflecting 
a different preferential transverse emission in the direc- 
tion of out-of-plane or in-plane, respectively. Meanwhile, 
absolute flow values for photons are smaller than the pro- 
ton's like the behavior of transverse flow. Except of the 
opposite sign, we sec that both elliptic flows have similar 
tendency with the increasing of px, i.e., their absolute 
values increase at lower px, and becomes gradually sat- 
urated, especially for direct hard photons. 

To explain the correlation of the collective flow between 
direct hard photons and free protons indicating above, we 
should note that flows of direct hard photons originate 
from the individual proton-neutron collisions. As Eq. [3] 
shows, we can roughly consider that in the individual 
proton-neutron center of mass system, in directions per- 
pendicular to incident proton velocity, i.e. 9-y = 90°, the 
possibility of hard photon production is much larger than 
the parallel direction, actually it is also in agreement with 
the theoretical calculations and the experiments [40, |41| . 
As a whole, the flow of hard photons should be corre- 
lated with the collective movement of the nucleons, and 
presents an opposite behavior. Consequently, flows of 
hard photons and of free nucleons exhibit an opposite 
character. 



IV. DEPENDENCES OF HARD PHOTON 

PRODUCTION AND FLOW ON SOME 

VARIABLES 

A. Impact parameter dependence 

It is well known that anisotropic flow mainly origi- 
nates from the initial asymmetric overlap zone of collid- 
ing nuclei which induces different pressures or rotational 
collective motion of the participant region and leads to 
anisotropic emission of the particles. Peripheral collisions 
have more initial asymmetry in overlap zone than central 



> 

CD 



(a) direct photon 



r (b) free proton 



^ 



^^ 






(c) direct photon 



(d) free proton 



r^^ 



1 



-1 

-2 

< 

-3 M 

8 ^ 

6 

4 

2 





0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

b/bmax 

FIG. 8: (a) and (b) are the directed transverse flow parameter 
F of direct photons and free protons, respectively, as a func- 
tion of the reduced impact parameter for the reaction ''"Ca + 
Ca collisions at 60A MeV. (c) and (d) same as above but 
for elliptic asymmetry coefficient V2- 
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FIG. 9: (a) and (b) is the directed transverse flow param- 
eter F of direct photons and free protons, respectively, as a 
function of beam energy for the reaction which is the same as 
Fig. [H but only semi-central events (40 — 60%). (c) and (d) 
same as above but for elliptic asymmetry coefficient «2. 



collisions, thus more anisotropic emission of the particles 
is expected. Therefore anisotropic emission shall be sen- 
sitive to impact parameter. 

Here we simulate the reaction ''^Ca + ''^Ca collisions 
at 60A MeV, and use the EOS parameters with the com- 
pressibility K of 235 McV for the nuclear mean field [/. 
Fig. m (a) and (b) shows the directed transverse flow pa- 
rameter F of direct photons and free protons, respec- 
tively, as a function of the reduced impact parameter (i.e. 
normalized by the maximum impact parameter) . Fig. [8] 
(c) and (d) are the same as (a) and (b) but for elliptic 
asymmetry coefficient vi- We can see that both F and v-i 
of direct photons and free protons have the similar ten- 
dency with impact parameter, i.e. their absolute values 
increase with the impact parameter except of slightly de- 
creasing in very peripheral collisions for V2. We also see 
that at all impact parameters, in contrast to to the neg- 
ative directed transverse flow and positive elliptic flow of 
free protons, direct photons show the positive F and the 
negative W2, i.e. the anisotropy is shifted by a phase 7r/2. 
It agrees with the previous conclusion that the azimuthal 
asymmetry of direct photons is anti-correlated with the 
corresponding free proton's flow. 



B. Beam energy dependence 

In Fig. [S] we show the directed transverse flow parame- 
ter F and elliptic asymmetry coefficient V2 of direct pho- 
tons and free protons, respectively, as a function of beam 
energy for the Ca -I- Ca reaction at (40 — 60%) central- 
ity. In the beam energy range studied here, same as the 
impact parameter dependence, the opposite signs of F 
and vi between direct photons and free protons are also 
anticorrelated. Moreover, except for the opposite sign, 
directed transverse flow parameter F of direct photons 
and free protons have similar structures with the increas- 



ing of beam energy. The value of direct photon ellip- 
tic asymmetry coefficient v-2, increases with beam energy 
from negative to positive, a tendency similar to avail- 
able experimental results on hard photons in relativistic 
heavy-ion collisions. 



C. EOS and symmetry energy dependences 

In this section we will discuss the nuclear equation 
of state and its symmetry energy term depdcnces of 
hard photon production. In our calculations, the single- 
particle potential taken as the Skyrme parameterized 
mean field potential including symmetry potential is 
shown below 

c/"^''^(p) = «(-)+/3(-r + K"lfHp,<5), (11) 

Po Po " 

where the coefficients a, h and a are parameters for nu- 
clear equation of state which is determined by the nu- 
clear saturation property and the incompressibility K of 
symmetric nuclear matter; 8 is the isospin-asymmetry 
parameter, b = {pn - Pp)/{pn + Pp), and Vasy {p,S) is 
the symmetry potential of neutrons (protons). 

First, we only investigate the effects of the incompress- 
ibility K on the hard photon production in symmetric 
reaction system. In the Skyrme potential, we take the in- 
compressibility K as 200 MeV, 235 McV and 380 McV, 
in which the first two correspond to soft and semi-soft 
potential and the last one corresponds to hard potential 
as we introduced in Sec. II A. 

Because of the sensitivity to the density oscillations of 
colliding system, hard photon should be also dependent 
of the EOS of nuclear matter, especially for thermal hard 
photons [a, 0] . Actually as shown in Fig. [TUl photon pro- 
duction shows its sensitivity to the compressibility, espe- 
cially for thermal photon which arc produced after t ^ 80 
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FIG. 10: (Color online) Time evolution of production rate of 
breinsstrahlung photons with time for inclusive events of *''Ca + 
'^"Csi at 30A MeV. Different EOS are used, namely the hard EOS, 
the semi-soft EOS and the soft EOS, respectively. 



FIG. 11: (Color online) Time evolution of hard photon mul- 
tiplicity with different incompressibility K for the reaction 
■'OCa + «Ca at 60A MeV. 



fm/c. It shows that the stiffer the EOS, the higher mul- 
tiphcity the thermal photons. In contrast, dnect pho- 
tons are produced by the first channel neutron-proton 
bremsstrahlung, their production rate only weakly de- 
pends on EOS since direct hard photons are emitted 
by the first channel n-p scattering when the system is 
in a highly nonequilibrium state during the compression 
stage, and they don not have enough times to feel the 
EOS. 

Shown in Fig. [Tl]is the time evolution of inclusive hard 
photon multiplicity from the symmetric reaction system 
"^"Ca + ''°Ca at 60^ MeV. We can see that at this en- 
ergy, for different incompressibility K, multiplicities of 
hard photons produced in the early stage of collisions are 
nearly equivalent and later they show a clear correlation 
with incompressibility K. Corresponding to direct and 
thermal photons, this indicates that direct photon pro- 
duction is not sensitive to incompressibility K, because 
at this reaction energy, comparable with two-body inter- 
actions, mean field can be neglected in producing hard 
photons in the early stage of collisions. However, there is 
a correlation between K and thermal photon production, 
with larger K producing more thermal photons. 

For the magnitude of flow parameters F and V2 , com- 
parisons are made for different equation of state. Fig. [12] 
shows the F and V2 as a function of impact parameter 
with different EOS. Generally, the directed flow becomes 
larger with the increasing of stiffness of the EOS, whereas 
the elliptic flow becomes smaller. 

In order to study the effects of symmetry potential, 
we now set incompressibility K as 380 MeV. We mainly 
consider three kinds of symmetry potential: (1) the sym- 
metry potential is neglected, that is 

V^°l{p,6)^0- (12) 

(2) the symmetry potential is a linear function of isospin- 
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FIG. 12: (Color online) Directed flow parameter (F) and elliptic 
flow {v2) as a function of impact parameter with different EOS 
parameters for Ca -|- Ca at 30A MeV. 



asymmetry parameter S, that is 



V2lip,S)^CsyraST-^P\ (13) 



where Csym=32 MeV, t"=1, tP=-1; 
(3) the Li et al. single-particle potential from the 
parametrization of nuclear symmetry energy used in 
Ref. [421 fo'^ studying the properties of neutron stars, 
that is Esymip) — Eo{po)i — )'^ ^ Li et al. derived a single- 
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FIG. 13: (Color online) Different symmetry potential which is 
used in the present model calculation as a function of density for 
an isospin asymmetry of 5 = 0.2 and the 7 parameters of 0.5 (red 
line) and 2.0 (green line), respectively. Solid lines are for neutrons 
and dashed lines for protons. 



FIG. 14: (Color online) Time evolution of hard photon mul- 
tiplicity with different symmetry potentials for the reaction 
■**Ca + **Ca at 60 A MeV. 



particle symmetry potential as [43|: 



V^?y{p.S) 



Eo{po){l- 



Po 



l)(f) 

p 



.4.2(-^)2/3 
Po 



12.7(^)2/3 

Po 



Sr"-ip), 



(14) 

where r"=l, tP—-1, Eo{po) sets as 30 MeV in this paper, 
and parameter 7 represents the stiffness of the symmetry 
energy. In the following, we consider two cases of 7 = 0.5 
and 2, respectively, corresponding to soft and stiff sym- 
metry energy, to explore the effects of different symmetry 
potential. 

The symmetry potentials for neutron and proton with 
S ~ 0.20 is plotted in Fig. [131 where the red line repre- 
sents Vasy, and the blue and green lines represent Vasy 
for 7 parameters of 0.5 and 2.0 for neutrons (upper 
branch) and protons (lower branch), respectively. For 
small isospin asymmetry and density near po the above 
symmetry potentials reduce to the well known Lane po- 
tential which varies linearly with i5 [ij] . Generally, the re- 
pulsive/attractive symmetry potential for neutrons/ pro- 
tons increases with density. 

Fig. [T4l shows the time evolution of inclusive hard pho- 
ton multipHcity in the reaction ''^Ca -I- ^^Ca at 60 A MeV. 
In order to separate the effects of symmetry potential 
from the incompressibility, we set K as a constant 380 
MeV. Fig. [14] presents the effects of different symmetry 
potentials on hard photon productions. Same as the in- 
compressibility K, direct photon production is insensi- 
tive to symmetry potential at the studied reaction en- 
ergy. For thermal photon, the curves of Vasy and V}sy 
have the nearly same multiplicities with time evolution. 



(2)/ 



.(2), 



and they are intermediate between Vasy {0.5) and Vasy (2). 
We note that thermal photon production is sensitive to 7- 
parameter in symmetry potential Vasyi'y), i.e. less 7, cor- 
responds to softer symmetry energy, induces more ther- 



mal photons. And we also note that the sensitivity of 
thermal photon production to 7 is not remarkable in 
comparison with its sensitivity to the incompressibility 
K, which is of course understandable since the symme- 
try energy is relative small term in comparison with the 
incompressibility K. 

From the effects of the incompressibility and symmetry 
potential on hard photon production, it appears that the 
yield of thermal photons has rather strong dependence 
on nuclear compressibility as well as symmetry energy, 
whereas the yield of direct photons is rather insensitive 
to nuclear EOS. Therefore thermal hard photon could 
serve as a porble of nuclear EOS in intermediate energy 
HIC. 



V. PHOTON-PHOTON CORRELATIONS 

In above sections, we mainly focus on the properties of 
inclusive hard photons, from which we find that produc- 
tion distribution and anistropic flows of hard photons 
reveal rich information of nuclear dynamics. Also, the 
sensitivity to the nuclear equation of state and symme- 
try energy for hard photons has been discussed. In this 
section, we will discuss the correlation properties of pho- 
tons in terms of two-particle correlation technique that 
can provide us a very powerful tool to characterise the 
properties of a particle source. In particular, two-boson 
relative momentum distributions enable one to determine 
the space-time structure according to the formalism of 
Bose-Einstein correlations. The magnitude of the corre- 
lation can be related to ther space-time distribution of 
the boson source. In the foUwoing calculations, we will 
construct two-photon momentum correlation function as 
well as azimuthal correlation function, from which the 
space-time structure of the photon souce and anistropy 
property can be indicated. 



A. Two-photon momentum correlation 

Intensity interferometry (also called HBT correlation) 
is used as a universal tool to study the properties of any 
boson sources such as stars [43|, or photon and meson 
sources in the early phase of heavy-ion collisions [4y| . The 
formalism was developed starting from optics and quan- 
tum statistics and was finally adapted to the dynamics of 
heavy- ion collisions |47H55| . We have performed the cal- 
culation which evaluates the correlation function directly 
from the photon source distribution given by the BUU 
calculation. We store for each ith pn collision its posi- 
tion ri and the associated photon probability distribu- 
tion Pi{ki). After the completion of the calculation, we 
analyze this output data to construct plane waves with 
four-momentum ki at r,; and calculate the two-photon 
probability for z ^ j as [3^, [5^] : 



Pl2 =Pi 



^gi(/si-?t+fc2-?l) I J^gi(ki-r]+k2-rt) 



::7> ^\ 



Pi(g)2[l + 2ABcos[g(?r-r 



(15) 
(16) 



where Pi 2 represents the probability to produce a 
pair without correlation, and A and B are the ampli- 
tudes related to the normalized probabilities of the direct 



P.(fcl)P,(fc2) 



and cross terms 



,-^ 



Pjih)P,{k2) 



This 



corresponds to Fourier transforming the photon source 
event by event. We set the weight of the interference 
in Eq. [16] as k, that is k = 2AB, which was predicted 
between 0.5 and 1.0 depending on the anisotropy of 
hard photon emission. The exact experimental filter 
was finally applied to the projection onto the Lorentz- 
invariant relative four momentum Q = -^/q^ — q^ of the 
resulting distribution P12 and Pi 2 of Eq. 1161 then 
the two-photon correlation function was calculated as 

CMQ) 



l + Accos[g(?T-F;)]==/(Q) 



For convenience, k was set to 0.75 for our calculations 
in order to take into account the established anisotropic 
component in the angular distribution of hard photons. 
The two-photon correlation function has been fitted with 
a Gaussian parameterization with correlation strength A 
and radius parameter Rq: 



C12 = 1 + \exp{-Q^R^Q) 



(17) 



In the above relation Rq is the space-time parameter 
conjugate to Q, which measures an invariant length de- 
pending on the source-size parameters R and r: 



Rq 




jr/Rnqo/qy 
1 - {qo/qr 



(18) 



One can easily know that, since qo ^ q, Rq is a first 
order measure of the spatial extent of the source, that is 
R « Rq . We can then calculate from Rq the rms radius 
of the source as the one of a static three-dimensional 
Gaussian source: Rrms — vSRq. 
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FIG. 15: (Color online) Two-photon HBT correlation function for 
the reaction *^Kr + ^*Ni collisions at 60^4 MeV in the laboratory 
frame, Ej > 25 MeV. The green squares represent for direct pho- 
tons, blue triangles for thermal photons, and red open circles for 
inclusive photons. The black closed circles represent for the exper- 
imental data taken from Ref. 1381. l39ll. 
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FIG. 16: (Color online) Two-photon HBT correlation function for 
the reaction ^''Ca -|- ^''Ca collisions at 60 A MeV. The symbols are 
the same as Fig. 115! and the dashed lines are fitting functions as 
Eq.[ni 



For the comparisons with experimental data, we calcu- 
lated the reaction ^^Kr -I- ^^Ni at 60^ MeV, employing 
here the following filter which is very similar to the exper- 
iment: Ej{l,2) > 25 MeV, detector positions between 
polar angles of 35° and 165° (orienting downstream), az- 
imuthal opening angles of 0±28° and 180 ±28°, and 18° 
for the minimum opening angle. Fig. [15] presents HBT 
correlation functions of direct photons (green squres), 
thermal photons (blue triangles) and inclusive photons 
(red open circles), respectively. We find that the cor- 
relation function of direct photons is much larger than 
thermal one, and the correlation function of inclusive 
photons is intermediate. This result is reasonable, be- 
cause direct photons are emitted at the early compressed 
stage of collisions, so they have stronger interference than 
thermal photons which are produced in later thermaliz- 
ing stage. We also see that our result can well reproduce 
the experimental data [33 JS^ , especially for the corre- 
lation function of thermal photons, which agrees with 
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the experimental function even the oscillation structure. 
Therefore, we can successfully reproduce the HBT corre- 
lation function of hard photons by the BUU simulation 

In order to further study the HBT correlation of hard 
photon and extract photon source information, we calcu- 
lated the symmetric reaction system ^°Ca -I- ''"Ca at 60A 
MeV in the CM. frame, and only the events of central 
collision (0 < bred < 0.2) were taken for simplification. 
In Fig. [ini we find that the correlation function of direct 
photons is the largest, the thermal one is the least, and 
the correlation function of inclusive photons are between 
them. Moreover, in the conditions of above reaction, all 
of them can be well fitted by Eq. [171 After fitting cor- 
relation function, we can obtain two useful fitting pa- 
rameters: correlation strength A and radius parameter 
Rq. As the results on top right corner show, direct pho- 
tons have the largest correlation strength A, the second 
is inclusive photons and the least one is thermal photons. 
Thus A is a parameter which is sensitive to the intensity 
of interference. And we also get three radius parameters 
Rq. Rq of thermal photons is nearly twice than direct 
photons', and inclusive photon is equivalent to the later. 
To explain this, we know that direct photon is mostly 
emitted in the early stage when the reaction system is 
strongly compressed, so the emission source of photons 
is small. And later on, the thermalizing system extends 
very much in company with the production of thermal 
photons. In this reaction energy at Eiab = 60A MeV, 
most of hard photons are produced in the early stage 
of the collision, so the spatial source extent of inclusive 
photons should be approximate to the direct photons. 
As a result, they have equivalent radius parameter Rq. 
Then we can get the rms radius of the photon source by 
the equation Rrms = v5Rq. Therefore, two-photon cor- 
relation function provides the information of hard pho- 
ton source, which is available to investigate the emission 
source during the collisions. 



B. Tw^o-photon azimuthal correlation 

From individual-photon azimuthal distribution, we 
find that direct hard photons exhibit the azimuthal asym- 
metric emission in intermediate energy heavy-ion colli- 
sions, especially negative elliptic flow parameter V2. Ac- 
tually it is difficult to extract the elliptic flow parame- 
ter V2 by the method of reconstructing reaction plane in 
the experiment, so as below, we will discuss to use the 
method of two-photon azimuthal correlation to extract 
the elliptic flow of direct photons directly. 

For particles in the same class, we defined the par- 
ticle azimuthal correlation function following the multi- 
fragment azimuthal correlation method [2^, [53, [S^ . The 
multi-fragment azimuthal correlation function is defined 
as follows: 



CiA^) 



Ncor{A(t>) 
NuncoriAct)) 



(19) 
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FIG. 17: (Color online) Top panel: Two-photon azimuthal cor- 
relation function for the reaction *'^Ca + ^'^Ca at 60^4 MeV for 
semi-central events (40%-60%), and the dashed line is the fitting 
function as Eq. 1201 Bottom panel: Individual-photon azimuthal 
distribution for the same reaction, and the dashed line is a fit of 
Fourier expansion. 



where Ncor is the A^ distribution for fragment pairs from 
the same event and Nuncor is the At/) distribution by ran- 
domly selecting each member of a fragment pair to form 
mixed events. The Ac/) between all selected fragments in 
an event are used to construct the correlation function, 
that is, "^2 '^'^ angles for n fragments. 

In our calculations for two-photon azimuthal corre- 
lations, we get Ncor from the A(/) distribution of pho- 
ton pairs in the same event and Nuncor from the mixed 
events. For these correlations, we may make a fit of the 
Fourier series with the expression 



C{A<P) = A[l + Ai cos(A(/)) + As cos(2A(/))], 



(20) 



where Ai and A2 arc treated as fit parameters. 

Under the assumption of statistically independent 
emission of particles with the same azimuthal distribu- 
tion F{(f>) in an event, the azimuthal correlation function 
is simply related to F{(j)) via the convolution 



C(A</)) 



27r 



F(</))F((/) + A(/))d(/). 



(21) 



One the other hand. F{(f>) can be described by the 
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FIG. 18: Excitation function of V2 for reaction '"'Ca + '"'Ca at 
60^4 MeV with centrality 40%-60%, open circles represent the am- 
plitude of V2 extracted from two-photon azimuthal correlation and 
red ones from the fit of Fourier expansion to individual-photon az- 
imuthal distribution. 



Fourier expansion: 



F{cj)) = ^ ocl 



2^u„cos(n(/)), 

n=l 



(22) 



where the different nth harmonic Fourier expansion coef- 
ficient Vn is defined as the nth anisotropic fiow, of which 
vi is directed flow parameter and V2 is so cahed elhptic 
flow parameter. 

If we only take vi and V2 terms in Eg. 1221 and substitute 
Eq. [22] into Eq. [211 we can derive the form of C(A^) 
as foflows: 

C{A(t)) = B[l + 2v\ cos{A(t>) + 2vl cos(2A</>)] (23) 

From Eq[2I] and Eal^ we then get the relationship 
between the fltted parameter A2 and the elliptic flow pa- 
rameter Vi |29| : 



«2 



VA272. 



(24) 



Now we make the comparison of elliptic flow parameter 
V2 extracted from two-photon azimuthal correlation with 
fitting value from individual-photon azimuthal distribu- 
tion. In Fig.[T7l we calculated the reaction ■*°Ca + '"^Ca 
at 60^ MeV for semi-central events (40%~60%) with the 
compressibility K = 235 MeV, and we extracted the ellip- 
tic flow parameters V2 by the two distinct methods above. 
Then we can sec that the absolute values of f 2 are nearly 
equivalent, except that we can only get the amplitude of 
elliptic flow from two-photon azimuthal correlation, and 
the fltting value of V2 from individual-photon azimuthal 
distribution is negative. Further, as shown in Fig. llSi we 
get the excitation function of 1)2 by the above two meth- 
ods which shows nearly same amplitude of V2 in various 
reaction energy but different signs. Therefore, in view 



of the difficulty to reconstruct reaction plane in exper- 
iment, two-photon azimuthal correlation provides us an 
alternative method to extract the amplitude of elliptic 
fiow parameter V2 of hard photons. 



VI. SUMMARY 

In summary, we have systematically investigated hard 
photon production by the process of proton-neutron 
bremsstrahlung as well as the behavior of azimuthal 
asymmetry of hard photons and free protons produced 
in intermediate energy heavy-ion collisions in the frame- 
work of the BUU model. As discussed in previous studies, 
hard photons can be separated into direct hard photon 
which is produced by the first channel neutron-proton 
bremsstrahlung and thermal hard photon which is pro- 
duced by the later stage neutron-proton bremsstrahlung 
during the system's evolution towards thermalization. 
The kinetic energy spectra of two kinds of photons show 
the exponential form which provides us with information 
on the apparent temperature in different stage of heavy- 
ion collision. 

The azimuthal asymmetry parameters have been in- 
vestigated in detail. The time evolution of directed flow 
and elliptic flow of hard photons exhibits rich structure 
as the density of the system oscillates during the pre- 
equilibrium and thermalization stage of reaction system. 
This structure indicates that the azimuthal asymmetry 
evolves with time. A nonzero directed transverse flow 
and elliptic fiow parameter have been predicted for the 
direct hard photons produced by the first channel proton- 
neutron bremsstrahlung process in the intermediate en- 
ergy heavy-ion collisions. The asymmetry parameters of 
hard photons are plotted as a function of rapidity and 
transverse momentum, which show contrary signs in com- 
parison with the flow of free protons, i.e. the azimuthal 
asymmetry of direct hard photons seems to be anticorre- 
latcd to that of the corresponding free proton's. There- 
fore we can expect the direct hard photon can server as 
a good probe to nuclear matter properties in the early 
stage of HIC. 

Different variables dependences are investigated for the 
hard photon production and anisotropic flow. We can see 
that the absolute value of both directed fiow and ellip- 
tic fiow of direct photons and free protons increase with 
the impact parameter, indicating that the fiow mainly 
refiects the initial geometric asymmetry of the collision 
zone in a given beam energy and revealing a rich struc- 
ture of fiows with increasing beam energy. For directed 
fiow, direct photons reach a minimum of around 40^ 
MeV where the free protons approach the maximum neg- 
ative value, while for elliptic flow, the absolute values of 
direct photons show decreasing trends with increasing 
beam energy. The EOS dependence for hard photons 
at different time indicates that the direct photon is not 
sensitive to the nuclear incompressibility nor to the sym- 
metry energy. However, for the thermal photon, its mul- 
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tiplicity increases with nuclear incompressibility. For a 
given nuclear incompressibility, the soft symmetry energy 
favors thermal photon production. 

Finally, we calculated the two-photon correlations, in- 
cluding HBT momentum correlations and azimuthal cor- 
relations. From two-photon HBT correlations, we can 
extract photon source information, such as intensity in- 
terference and the spatial extent of the emission source. 
We also find that two-photon azimuthal correlations can 
provide us with a good method for extracting the ampli- 
tude of the elliptic flow parameter V2 of hard photons in 
the experiment. 

In light of the present study, we expect that direct pho- 
tons would be a very useful probe for exploring nuclear 



reaction dynamics in intermediate energy heavy-ion col- 
lisions, while the thermal hard photon can give us some 
hints on the nuclear EOS, including the symmetry en- 
ergy. 
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